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Abstract 

Cosmic rays represent one of the most fascinating research themes 
in modern astronomy and physics. After almost a century since their 
discovery, a huge amount of scientific literature has been written on 
this topic and it is not always easy to extract from it the necessary 
information for somebody who approaches the subject for the first 
time. This has been the main motivation for preparing this article, 
which is a concise and self-contained review for whoever is interested 
in studying cosmic rays. The priority has been given here to well 
established facts, which are not at risk to get obsolete in a few years 
due to the fast progress of the research in this field. Also many data 
are presented, which are useful to characterize the doses of ionizing 
radiation delivered to organisms living on the Earth due to cosmic 
rays. The technical terms which are often encountered in the scientific 
literature are explained in a separate appendix. 



1 Introduction 

Cosmic rays (CR) represent a fascinating subject of research, which is re- 
cently witnessing a growing interest within the scientific community. With 
the next generation of cosmic ray detectors, like for example the Pierre Auger 
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Observatory which is currently under construction, there is the hope that the 
questions posed by these cosmic particles will find soon some satisfactory an- 
swer. One should also note that the attention to CR is not only restricted to 
the traditional fields of high energy physics and astroparticles. This article 
is for instance the result of the authors' efforts to investigate the mutagenic 
effects of cosmic rays on cells. To this purpose, one needs to characterize 
very precisely the fluxes and intensities of particles arriving to the ground as 
an effect of the interaction of CR with the atmosphere of the Earth. 

The increasing popularity of CR is accompanied by an increasing demand 
of information on this topic, but it is not always easy to extract this informa- 
tion from the huge amount of literature which has been written on CR during 
a period of almost a century after their discovery. An additional problem is 
that, looking at the scientific literature on CR, often concepts like the integral 
vertical intensity or the differential integrated flux are encountered. These 
terms sound somewhat puzzling and exotic for a reader that encounters them 
for the first time. These motivations have compelled us to prepare this short 
review, which would like to be a compact but self-contained reference on CR. 
The preference has been given here to well established facts, which are not 
in danger to become obsolete in a few years due to the fast progress in this 
topic. Moreover, a considerable effort has been done in order to explain in 
details and with the help of figures the technical terms used in the current 
scientific literature. The puzzles raised by the existence of ultra-high energy 
cosmic rays with energy of 10 19 eV or higher, like the mystery of their ori- 
gins or the apparent violation of the theory of relativity which is connected 
with them, have been just briefly mentioned. Hopefully, these puzzles will 
be solved with the next generation of cosmic rays detectors. Another aim of 
this work is to present data concerning the physical parameters, e. g. types 
of radiation, delivered effective doses and dose rates, fluxes and intensities of 
incoming particles, which characterize the radiation to which the organisms 
living on the Earth are exposed because of CR. These data are certainly 
of interest for scientists working in life sciences. Cosmic rays are in fact the 
source of an almost uniform background of ionizing radiation which is present 
everywhere on the Earth. Most of their energy arrives to the ground in the 
form of kinetic energy of muons. The latter particles are very penetrating 
and are able to travel for kilometers in water and for hundreds of meters in 
rock. Since ionizing radiation is mutagenic, it is very likely that the radiation 
of cosmic origin has shaped in some way the evolution of life on our planet, 
generating some kind of adaptive response in cells. 
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Finally, as already mentioned, this article is self-contained, but of course 
it is far from being complete, because it is impossible to cover all the immense 
literature which exists on cosmic rays. To integrate the material presented 
here, suggested further readings are for instance [HElinillllHllEl and references 
therein. 

2 Cosmic rays and the natural background ra- 
diation 

2.1 Introduction to cosmic rays 

Cosmic rays, first discovered by Victor Hess in 1912, are charged particles 
accelerated at very high energies by astrophysical sources located anywhere 
beyond the atmosphere of the Earth. 89% of cosmic rays consists of protons, 
followed by a— particles (~ 10%) and heavier nuclei (~ 1%) 1 . All elements 
of the periodic table, including transuranic elements, have been detected. 
More details on the composition of cosmic rays can be found in Ref. [2] and 
references therein. 

Within the flux of incoming particles one can distinguish different com- 
ponents. The most relevant ones are galactic CR (GCR), solar energetic par- 
ticles (SCR) and the anomalous CR (ACR) |Hj. GCR originate from sources 
located in our galaxy outside the solar system. It is believed that GCR are 
a consequence of astrophysical events like stellar flares, stellar coronal mass 
ejections, supernova explosions, particle acceleration by pulsars [Zj. Very of- 
ten, the word cosmic rays refers only to GCR [T]. The smallest detectable 
energy of GCR is about 1 GeV. Below this energy, the screening effect of 
the solar wind is too strong to allow them to penetrate the heliosphere p]. 
A detailed information on the interaction of GCR with the magnetic fields 
of the Sun and of the Earth is contained in Ref. [2j and references therein. 
Here we would just like to mention the action on CR of the galactic magnetic 
fields which are generated by the spinning of the Milky Way. These fields 
are relatively weak, because the average magnetic field in our galaxy is of 

1 These data are taken from Ref. Of course, as it always happens in the case of 
experimental data, there is some uncertainty due to measurement errors. For this reason, 
different authors report slightly different values as in the case of [Jj, which gives 95% for 
protons, 3.5% for a particles and 1.5% for all the rest. It should also be noted that the 
flux of particles below the energy of lOGeV strongly depends on the 11-year solar cycle. 
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the order of 10 T. However, since CR consist of charged particles traveling 
along huge distances, at the end these magnetic fields are able to bend their 
trajectories in a relevant way. At this point one should note that there are 
two different components in the magnetic fields of our galaxy, a regular one 
and a turbulent one, see Ref. [TH]. The strengths and directions of the mag- 
netic fields belonging to the turbulent component are random [HI EH- Due 
to this randomness, the trajectories of CR are randomly bent. As a conse- 
quence, the flux of CR which arrives on the Earth is also random or, more 
precisely, isotropic. For this reason, it is not easy to ascertain where GCR 
are coming from. Besides being isotropic, the particle flux is constant in time 
too, so that GCR form an almost uniform background of ionizing radiation 
on the surface of the Earth. CR of energies up to 10 21 eV have been observed. 
These are considerable energies for a microscopic particle. For example, the 
upper energy limit of 10 21 eV corresponds in SI units approximately to 160 
Joules. This is comparable to the kinetic energy of a ball of 0.8kg thrown at 
the speed of 50km/hour. The origin of such ultra high energy CR (UHECR) 
is so far unknown, but there are strong hints that they could be produced 
outside of our galaxy. Candidate sources of UHECR could be relativistic 
plasma jets from supermassive black holes [TT], explosions of galactic nuclei 
[Zj, but other possibilities have been proposed, like magnetic monopoles, see 
for example [T2]. The main evidence which suggests that UHECR are of ex- 
tragalactic origin is provided by the fact that the magnetic fields which are 
present in the Milky Way are not able to trap CR of that energy. Indeed, 
already protons of energy higher than 10 15 eV are able to escape the galactic 
confinement 2 . As a consequence, if protons of energies of 10 19 eV or higher 
would be produced by sources located in our galaxy, they would escape it in 
all possible directions following trajectories which are almost straight lines. 
For this reason, the ultra high energy protons which reach the Earth should 
arrive along directions which are approximately parallel to the galactic plane. 
However, this conclusion is not confirmed by observations, see Refs. [THl EH 
and references therein. Observations show in fact that the spatial distribu- 
tion of ultra high energy protons is isotropic, so that their directions are not 
aligned with the galactic plane. This is of course a strong hint that UHECR 
are of extragalactic origin. 

2 In the case of heavier nuclei, the threshold energy for escaping the galactic confinement 
is higher than that of protons. It is for this reason that one observes a greater proportion 
of heavier nuclei with respect to protons in CR with energy above 10 15 eV. 
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On the other side, there is at least one argument which seems to point 
out that the sources of UHECR are not very far from our galaxy. In fact, 
it has been noted that protons of energies above 5-10 19 eV would lose their 
energy by interacting with the photons of the microwave (big bang) back- 
ground. This effect was predicted in 1966, one year after the detection of the 
microwave background radiation, by Kenneth Greisen, Vadem Kuzmin and 
Georgi Zatsepin. The energy threshold of 5-10 19 eV is called the GZK-limit, 
from the names of its discoverers. Protons with energies above that threshold 
will be slowed down during their travel to the Earth by the mechanism of 
energy-loss pointed out by Greisen, Kuzmin and Zatsepin until their energy 
falls below the GZK-limit. This mechanism is so efficient that, in practice, 
protons with energies higher than 5-10 19 eV should not be observed on the 
Earth if their source is located at distances which are greater than 50Mpc 3 . 
Since ultra high energy protons have instead been detected, this implies that 
they originate from sources which are within the range of 50Mpc. Yet the 
known cosmic objects which could be able to accelerate protons to such high 
energies are at least at a distance of about 100 Mpc or more. So cosmic 
ray protons above this energy should not arrive on the Earth or we should 
explain why their formidable sources, which should be relatively near to us, 
remain invisible. These contradictions are known under the name of GZK 
paradox. On these points see Ref. [15J. 

Let's now discuss the remaining components of CR. Energetic solar events, 
like for instance solar flares, are able to accelerate particles up to some GeV 
very efficiently within the time of 10 seconds. The SCR are mainly pro- 
tons, heavier nuclei and electrons. Finally, it is worth mentioning also the 
case of ACR. This kind of CR is mainly characterized by ions of elements 
which are difficult to ionize, including He, N, O, Ne and Ar. Moreover, ACR 
have a relatively low energy, up to a few hundreds of MeV [16J. It is thus 
improbable that CR of such a low energy could originate from the violent 
phenomena which produce GCR. Besides, we recall that CR of energy be- 
low lGeV coming from outside the heliosphere cannot penetrate very deeply 
inside the solar system, because they are deflected by the solar magnetic 
fields. There are indeed evidences that ACR may be neutrally charged dust 
particles present in the interstellar gas near the border of the solar system 

3 Mpc stands for megaparsec. 50Mpc are approximately 150 millions of light years. This 
is about 1500 times the diameter of a galaxy and it is not a big distance in comparison 
with the cosmic scale of distances. 
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[H El El- When these particles enter in contact with the far edges of the 
heliosphere, they are ionized by UV solar photons or by interactions with 
the particles within the solar wind. Once these dust particles have been ion- 
ized, they are accelerated by the shock waves formed when the solar wind 
encounters the interstellar plasma. Particles escaping the shocks may diffuse 
toward the inner heliosphere and arrive on the Earth as ACR. Actually, there 
are many open questions on the mechanism of ACR acceleration. Hopefully, 
these questions will be answered in 2007. By that date, in fact, the Voyager 1 
should reach the region in which the shocks occur, which is thought to be 
somewhere between 75 and 100 AU from the Sun 4 . This will be the first 
time that an example of CR acceleration will be observed directly. 

2.2 Interaction of CR with the Earth's atmosphere 

When CR arrive near the Earth, they hit the nuclei of the atoms of the atmo- 
sphere, in particular nitrogen and oxygen, producing in this way secondary 
particles. The first interaction of the CR primary particle takes place in the 
top 10% of the atmosphere [21] • The most relevant reactions 5 , remembering 
that 90% or more of CR consists of protons, are: 

pp — > pnir + or pp — > pp%° (1) 

pn — > ppn~ or pn — > pmr° or pn — ► nnn + (2) 



In the above reactions all the secondary particles are hadrons, namely protons 
p, neutrons n and pions in all their charged states 7r ± , 7r°. Pions may in turn 
decay according to the following processes: 



where the \i s are muons, the 7's are photons and u^, are respectively 
muonic neutrinos and their anti-particles. The mean life of pions is 26ns 

4 At the time in which this article was finished, a distance up to about 90 AU has 
been explored by Voyager 1. Some of the results of the observations can be found in 

Refs. da muni. 

5 One should remember that the collisions of CR with the atmosphere gives raise also to 
less relevant reactions, which produce particles like kaons, 77 particles and even resonances. 



(3) 
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(4) 
(5) 
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for ^ and 10~ 16 seconds in the case of it . For this reason, charged pions 
may still collide with air atoms before decaying, but it is very unlikely that 
this happens in the case of neutral pions, which have a very short average 
life. Other secondary particles, like protons, neutrons and photons interact 
very frequently with the atoms of the atmosphere giving raise in this way to 
a cascade of less and less energetic secondary particles. At the end, these 
particles are stopped by the atmosphere or, if the energy of the primary 
particle was sufficiently high, they can reach the ground. 

The main mechanism of energy loss 6 of high energetic hadrons consists 
in the disintegration of the molecules of the atmosphere |22], see Fig[TJ This 
leads to the creation of new particles through nuclear interactions like those 
shown in Eqs. (JTJ) and (j2J). At lower energies, dissipative processes become 
predominant, in which the molecules of the atmosphere get either ionized or 
excited. The most relevant process of this kind in the case of heavy charged 
particles is the ionization of the molecules of the atmosphere. Lighter charged 
particles like electrons and positrons lose their energies not only by ionization, 
but also by bremsstrahlung. This consists in the radiative loss of energy of 
charged particles moving inside matter, when they are deflected by the elec- 
trostatic forces of the positive charged nuclei of the surrounding molecules. 
Photons and neutrons, the remaining relevant particles in the cascade, are 
examples of indirectly ionizing radiation. Their interaction mechanisms are 
more complicated than those of charged particles and will not be described 
here. The interested reader may find a more detailed account on the way in 
which radiation of different kinds interacts with matter in Ref. [23] . 

The total number of secondary particles N sec within the cascade grows 
rapidly, mainly sustained by the processes of bremsstrahlung and pair pro- 
duction due to electrons, positrons and photons. Hadrons like protons and, 
to a less extent, neutrons, are easily stopped by the atmosphere, so that they 
increase relevantly the number of particles by disintegrating the molecules of 
air only during the first stages of the formation of the cascade. The decays of 
pions given in Eqs. (J3| and (JHJ) produce muons and photons of considerable 
energies. The muons are very penetrating particles and do not interact very 
much with the air. They lose a small fraction of their energy before reaching 
the ground by ionizing the molecules of the atmosphere. The photons give 
raise instead to electron-positron pairs e + e~ . In turn, electrons and positrons 
create other electrons by ionization or other photons due to bremsstrahlung. 

6 Here and in the following energy means the kinetic energy of the particles. 
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In this way, while the cascade propagates inside the atmosphere, the number 
of its electrons, positrons and photons grows almost exponentially. The max- 
imum number of particles inside the cascade is attained when the average 
energy per electron reaches the threshold Et ~80MeV. When the energy of 
electrons in air falls below that threshold, ionization starts to prevail over 
bremsstrahlung as the main mechanism of energy loss of electrons in air and 
the process for increasing the number of particles described above ceases to 
be effective. 

If the energy of the primary particle is below 10 14 eV, essentially only the 
penetrating muons and neutrinos are able to arrive to the sea level, while 
the other particles in the cascade are absorbed at higher altitudes. Actually, 
neutrinos interact so rarely with matter that they could pass through a light 
year of water without undergoing any interaction. Thus, if one is concerned 
with the dose of ionizing radiation delivered by CR to the population, the 
contribution of neutrinos can simply be neglected. Muons are more dangerous 
for the health. They have a short mean life at rest (2.2ns), but since they 
travel at very high speeds, they manage to reach the surface of the Earth 
due to the relativistic dilatation of time. Part of these muons can still decay, 
giving raise to electrons e~ or positrons e + , mainly according to the process: 
^ = e ± + v e + zv 

When the energy of the primary particles is instead above 10 14 eV, the 
cascade of secondary particles arrives to the ground before being stopped 
by the atmosphere. In that case, the cascade is called with the name of 
air shower, see FigHJ To be precise, the effects of an air shower started 
by a primary particle of 10 14 eV are relevant up to altitudes comparable to 
that of Mount Everest OEHl. Only air showers generated by primaries of 
energy of about 10 15 eV or higher are able to reach also the typical altitudes 
of inhabited areas and arrive to the sea level. The frequency of these air 
showers is relatively high, because the total flux of primary particles with 
energy E > 10 15 eV is of about 100 particles per m 2 per year. Giant air 
showers produced by primaries of energies beyond 10 20 eV are much more 
rare: their total flux is of 1 particle per km 2 per century [T]. More data 
concerning fluxes of incoming primary particles and an explanation of how 
these data are measured using ground detectors, can be found in Refs. [6|l26]. 

In the air shower we distinguish a nucleonic component, a muon compo- 
nent and an electromagnetic component. The nucleonic component is gen- 
erated by high energetic protons and neutrons, which disintegrate the atoms 
of the air giving raise to other protons and neutrons. The fluxes of electrons, 
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Figure 1: This figure illustrates schematically how air showers generate from 
cosmic rays. The high energetic primary particle, usually a proton, whose 
trajectory has been denoted in black, starts to interact with the molecules of 
the upper atmosphere. In this way, secondary particles are produced, which 
give raise to other particles (tertiary, quaternary etc.) via other interactions 
with the atmosphere or via decay processes. The total flux of particles can be 
divided in a electromagnetic component (photons, electrons and anti-electrons 
or positrons), see the blue trajectories in the figure, in a muon component, see 
the cyan and red trajectories, and finally in a nucleonic component (mainly 
protons, neutrons, rarely pions) denoted in brown. At the sea level, the air 
shower has the form of a pancake, whose height is around two meters, while 
its radius usually around a few hundreds of meters, but may reach some 
tents of kilometers in the case of very energetic cosmic rays. The nucleonic 
component is usually confined in a narrow cone centered along the direction 
of the incoming primary particle. 
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positrons and photons started by the decay of the 7r°'s, together with the 
electrons and positrons coming from the decay of muons or from the ioniza- 
tion due to the hadrons, form the electromagnetic component. In the early 
times of CR research, electrons and positrons have been called the soft com- 
ponent, while muons coming from the decay of the charged pions have been 
called the hard component. These names originate from the fact that muons 
are very penetrating and thus they may be regarded as "hard" particles. For 
instance, at the energy of 1 GeV, the range 7 of muons is 2.45-10 5 g cm -2 . This 
implies that in water, which has a density of 1 g cm -3 at 4 C, muons run 
along an average distance of 2.45 km before being stopped completely. In 
standard rock, which has a density of 2.65 g cm -3 [2], this average distance 
reduces to about 900 meters. 

At the sea level, the air shower has approximately the form of a pancake 
with an height of 1-2 meters. Its extension in the other two directions, 
defined as the distance in which 90% of the total energy of the shower is 
contained, is given by the so-called Moliere radius. For example, in the case 
of an air shower started by a primary particle of an energy of 10 19 eV (10 
EeV) the Moliere radius has a length of about 70 meters. The real extension 
of the shower is much bigger and some of the muons may be detected up 
to a distance of a few kilometers from the core [22] . Usually the nucleonic 
component, which is composed by heavier particles than those of the muon 
and electromagnetic components, is less deflected from the direction of the 
incident primary particle by the interactions with the atmosphere and it is 
thus concentrated in a narrow cone inside the air shower. The center of the 
cone is roughly aligned with the direction of the original primary particle. 
The number of secondary particles which arrives on the ground with an 
air shower is huge. Considering particles whose energies are greater that 
200keV, an air shower generated by a lOEeV primary particle contains up to 
10 10 particles, mostly photons, electrons and positrons. Electrons outnumber 
positrons in the ratio 6 to 1. The maximum number of particles, i. e. the 
so-called point of shower maximum or simply shower maximum, is attained 
at an altitude of 2-3 km above the sea level. Many other data and diagrams 
describing the propagation of air showers in the atmosphere may be found 
in Ref. [ZZ|. 

7 The range is defined as the average depth of penetration of a charged particle into 
a material before it loses all its kinetic energy and stops. The concept of range has a 
meaning only in the case of charged particles whose energy is kinetic energy which is lost 
continuously along their paths due to ionization and bremsstrahlung processes (2^1 • 
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When air showers approach the ground, about 85% of their energy is 
concentrated in the electromagnetic component. The contribution of the 
muon and nucleonic components is thus much less relevant. The situation 
changes completely if we consider all CR, and not only those which have 
sufficient energy to give raise to an air shower. As we see in Fig. [2 muons 
are in fact responsible for about 85% of the total equivalent dose delivered 
by CR to the population at the sea level. As a consequence, globally it is 
the muon component the most significant from the energetic point of view 
and not the electromagnetic component. The reason of this fact is that 
primary particles with energy E > 10 15 , namely those which can produce air 
showers, form just a minimal fraction of the total amount of CR arriving on 
the Earth. For example, the total flux of particles with energy E > 10 12 eV 
(1 TeV), is of 1 particle per m 2 per second, i. e. a factor of 3-10 5 higher 
than the total flux of CR with energy E > 10 15 reported above. In other 
words, there is an overwhelming number of CR with energy lower than 10 15 eV 
which are not able to start an air shower, but may still generate energetic 
muons. Being very penetrating particles, these muons are not stopped easily 
by the atmosphere and reach the sea level altitude, where they represent the 
biggest source of ionizing radiation of cosmic origin. Other particles which 
deliver relevant doses of ionizing radiation to the population on the surface 
of the Earth are photons, electrons and neutrons. The percent contributions 
to the total equivalent dose of the various components of CR as a function 
of the altitude is given in Fig. [21 In that figure, which has been published 
in 1996, the curve concerning neutrons should be taken with some care, 
because the data on neutron fluxes in the atmosphere were still sparse at 
that time [Hj. Other data about energies and fluxes of particles due to CR 
will be given in the next Section. We note that protons and neutrons prevail 
at higher altitudes, but they are rapidly absorbed by the atmosphere and 
muons become dominant at lower altitudes. 

2.3 Intensities and fluxes of CR 

CR are the source of an avalanche of secondary particles which hit constantly 
the surface of the Earth. To determine the energy and number of these par- 
ticles, their directions of arrival and their distribution in time, quantities 
like the integral vertical intensity or the differential directional intensity are 
measured. The meaning of these quantities is explained in details in a sep- 
arate Appendix at the end of this article. In this section, we present some 
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Figure 2: Percent contribution of the various CR components to the total 
equivalent dose at different altitudes. This figure is based on data of Ref. 

experimental data which are useful in order to characterize the contribution 
of the muon, electromagnetic and nucleonic components to the background 
of ionizing radiations on the ground due to CR. 

The integral vertical intensity (IVI in the Appendix, see Eq. (|2Tjl ) of the 
muon component with energy above 1 GeV at sea level is approximately 
I^ rd (9 = 0) ~0.70-10~ 2 cm-V^r- 1 ^. The integral directional intensity of 
muons in the other directions, which are at an angle 6 with respect to the 
vertical direction, has the following behavior: /^ rd (6') oc l-%f d (0 = 0) cos 2 9. 
More complete phenomenological formulas for the angular distribution of CR 
intensities may be found in Refs. OEHIEHI- As Fig. El shows, muons arriving 
to the ground are very energetic. The most frequent muon energy is 500MeV, 
while the average muon energy is 4GeV. There is almost no protection from 
this source of radiation, since, as we have seen before, high energetic muons 
are able to penetrate thick layers of concrete and rocks. 

The integral vertical intensity for electrons plus positrons with energies 
greater than 10, 100 and 1000 MeV is very approximately given by[2]: 

lf v t po \9 = 0, E min > 10MeV) ~ 0.30 • lO^cm-V^ -1 (6) 
lf v t pos (9 = 0, E min > 100MeV) ~ 0.06 ■ lO^cnrW- 1 (7) 
lf v f po \e = 0, E min > 1000MeV) ~ 0.02 • lO^cm" W" 1 (8) 

Moreover, the total flux of electrons plus positrons amounts approximately to 
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Particle fluxes at the sea level 




Particle energy (MeV) 

Figure 3: Differential integrated ffux (see Eq. fpTTT) of the Appendix) of the 
different components of CR related radiation at the sea level. This figure has 
been created on the basis of an analogous figure appeared in Ref. [HOI- The 
original data are taken from Ref. |32] . 

30% of the total particle flux which is reaching the ground due to cosmic rays. 
Always according to j2J 3 the ratios of photons to electrons + positrons is ap- 
proximately 1.3 above 1 GeV and 1.7 below the critical energy E T ~ 80MeV 
mentioned above. The differential flux on the ground of the various com- 
ponents of radiation related to CR can be seen in Fig. El One may observe 
that different particles become predominant at different energies. In the low- 
est portion of the energy range, electrons are predominant. At energies of 
around 1 MeV electrons are taken over by fast neutrons which arise mainly 
due to the de-excitation of atmospheric nuclei following compound-nucleus 
reactions pi] 8 . Electrons become once again predominant in the energy 
range going from a few MeV up to some tens of MeV. Starting from energies 

8 Roughly speaking, in a compound nucleus reaction a neutron or a proton, but also 
an a— particle interacting with the nucleus of an atom creates a nucleus of higher atomic 
number which is metastable and decays after a short period of time. For example, a 
possible compound-nucleus reaction is: p + 63 Cu — > 64 Zn*. The de-excitation of the 
metastable nucleus produces others neutrons and protons, e. g.: 64 Zn* — > 63 Zn + n or 
64 Zn* — > 62 Cu + n+p. Compound-nucleus reactions become possible only if the energies 
of the incoming nucleons or alpha particles are such that the de Broglie wave lengths of 
these particles are comparable with the size of the hit nucleus. 
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approximately above 200MeV, the number of muons becomes overwhelm- 
ingly high in comparison to that of the other particles. In considering the 
above data, one should of course take into account the fact that, at low ener- 
gies, let's say below lOMeV, particle fluxes are strongly dependent on many 
factors, including weather conditions, so that there are big uncertainties in 
their measurement up to an order of magnitude [32] • Moreover, all the data 
presented so far in this Section refer to the sea level altitude. With increasing 
altitudes, the contribution to the particle flux given by protons and neutrons, 
the nucleonic component of CR, becomes more and more relevant (see Fig. EJ) 
and is predominant above atmospheric depths of approximately 500 g cm -2 9 . 

Altitude vs atmospheric depth 
45 I 1 1 1 1 1 1 1 1 1 1 

<> 
40<> 




100 200 300 400 500 600 700 800 900 1000 
Atmospheric depth (g/cm 2 ) 



Figure 4: Conversion diagram from atmospheric depth to altitude. The data 
corresponding to the dots are taken from Ref. [33] • 

To conclude this Section, we provide also some data regarding the particle 
flux outside the heliosphere. The total flux of CR in the galaxy is large, 
about 100000 particles m~ 2 s _1 [31]. Much lower is for example the integral 
directional flux of CR primaries with E > 2 ■ 10 15 eV is about $ idf = 75000 

9 The atmospheric depth is a quantity which is often used to measure the altitude. For 
convenience, we give in Fig. a diagram which is useful to make the conversion from 
atmospheric depth in g/cm 2 to altitude in kilometers. 
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particles km -2 sr -1 day -1 . This datum has been derived using a formula for 
the integral directional flux given in Ref. [25J, which is based on the results 
of measurements and simulations of incoming cosmic ray fluxes reported in 
Ref. |35j. The dependence of this flux on the direction of the incoming 
particles is minimal since, as mentioned before, the distribution of these 
particles is isotropic due to the presence of random magnetic fields in the 
Milky Way. Finally, the integral directional flux of particles with energies 
above 10 20 eV is = 1 particle km -2 sr -1 century -1 [36j. With such a 
small flux, the investigation of CR with energy beyond the GZK limit requires 
detectors which cover a very large area. Perhaps with the next generation of 
detectors it will be possible to solve the puzzles connected with UHECR. 

2.4 Dose of ionizing radiation from CR in present times 

In present times 10 the effective dose rate (EDR) delivered by CR to the 
human population varies from a minimum of 300/iSv year -1 to a maximum 
of 2000/xSv year -1 . This wide range depends on many factors, the most 
important one being the altitude. At sea level, the value which usually is 
given for the EDR is 270 /iSv year -1 , or equivalently 31 nSv h -1 . One 
should keep in mind however that this estimation is the result of population- 
weighing average. As a matter of fact, even if the altitude is fixed at the 
sea level, the EDR still changes with the latitude within a range of variation 
of approximately 10%. Strictly speaking, 270 /iSv year -1 is the dose rate 
received by the population living at a latitude which is near the 30° parallel. 
It turns out in fact from the distribution of the population on the Earth that 
this is the average latitude at which people are living. Besides, the above 
value of EDR takes into account only the contributions of muons and of the 
electromagnetic component. The nucleonic component, which at the sea level 
is essentially consisting of neutrons, gives to the average EDR at the sea level 
an additional contribution of 48 /iSv year -1 or, equivalently, 5.5 nSv h . The 
figures concerning neutrons should be taken with some care, since up to the 
time in which Ref. [H] was released, the available data on neutron fluxes were 
sparse. If one considers also the different altitudes in which the population 
lives, the population- weighted EDR is of 380 /iSv year -1 , corresponding to 
an average habitant of the planet Earth living approximately near the 30° 

10 The data concerning the present levels of radiation coming from CR are taken from 
the United Nations Report of the year 2000 [Hj. 
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parallel and at an altitude of 900 meters above the sea level. 

3 Conclusions 

In this review a short but thorough account has been provided on what it is 
known about cosmic rays, starting from their origin in space and arriving to 
the doses of ionizing radiation delivered by them to the human population. 
Only the sources of CR have not been discussed here, because this argument 
is outside the aims of this article. Much attention has been dedicated to 
fluxes and intensities of CR and of the particles arriving on the ground as an 
effect of the cascades initiated by CR in the atmosphere. These quantities are 
of interest for scientists working in different subjects. Apart from research in 
high energy physics and astronomy, the fluxes and intensities of particles of 
cosmic origin are also studied for radioprotection purposes [3 El and for their 
capability of causing potentially harmful failures in computers and electronic 
storage devices |34j . Fluxes of CR are also carefully measured due to their 
relevance to space explorations, see for example [37J. Finally, in Appendix 
A the various kinds of fluxes and intensities of CR and related particles 
which one encounters in research articles about CR have been defined and 
their meaning has been illustrated. Concrete expressions for these quantities 
have been given in terms of mass densities, velocity distributions and energy 
densities. Both relativistic and non-relativistic cases have been treated. Up 
to now, a systematic classification and explanation of these quantities such 
as that provided in this work was missing in the scientific literature on CR. 
The necessity of filling this gap justifies the length of this Appendix. 

A Appendix A: Definitions of intensity, flux 
and related quantities 

We have seen that, in order to characterize the intensity and flux of charged 
particles which arrive on the Earth due to CR there exist an entire zoo of 
observables. Their names and meanings may sound puzzling for somebody 
who is not acquainted with them. Moreover, the same observable is some- 
times called with different names by different authors, or, on the contrary, the 
same name is used to describe two slightly different observables in different 
contexts. Besides, it is not easy to find an explanation of these observables 
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in the scientific literature. Books on radiative transport contain often useful 
information, see for example [38J, however these books describe the inten- 
sity and flux of radiation emitted by an energy source. Here we are instead 
dealing with intensity and flux of particles arriving at a detector. For these 
reasons and also to make this article self-contained, in the following it will 
be made an effort to explain the meaning of the various quantities which are 
relevant in the physics of CR. 

A.l Differential directional intensity 

The differential directional intensity (DDI) Iddi is defined 11 in such a way 
that the quantity 

dNi = I ddi dSdQdEidt (9) 

represents the number of particles of a given kind incident upon the infinites- 
imal element of area dS during the time dt within the element of solid angle 
dVL perpendicular to dS and within the energy interval [Ei, Ei+dEi]. Here the 
index % = 1, 2, 3, . . . labels the different kinds of particles (electrons, protons, 
muons etc.) 

To compute explicitly the Iddi in terms of physical parameters like particle 
velocity and mass or energy density, let us consider a point P in the space, 
whose position with respect to a cartesian system of coordinates Oxyz is given 
by the radius vector r = (x,y,z). In the following, it will be convenient to 
define a second reference system with origin in P and spherical coordinates 
zu, 9, (j), with < 9 < 7r and < <fi < 2-ri, see Fig. El We introduce also 
the infinitesimal vector element of surface dS = dSn. The area and the 
orientation of dS are given by dS and by the unit vector n which is normal 
to dS respectively. The element of surface dS is centered around the point 
P. Now we wish to count the number of particles of a certain type, e. g. 
electrons, which hit per unit of time the surface dS and whose velocities Vj 
are oriented according to a certain direction, given for instance by the unit 
vector e R (9,(j)). In mathematical terms this last condition is expressed as 
follows: Vj = |vj|efl(6 l , 0). We will see that, in the case of the flux, the 
direction of Br(0, (j)) may be arbitrary. However, in the case of the intensity, 
the element of surface dS should be by definition perpendicular to the vector 

11 Here we follow Ref. [221; in which a very clear and precise definition of the related 
concept of directional intensity is presented 
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e R (9, (ft). In other words: 

n = e R (9,<j>) (10) 

If the particles are non-relativistic, one may express the DDI in terms of the 
velocity and mass density of particles: 

v i = \v i (E i )\e R (e,<t>) (11) 

Pi = Pi (T,8,<f>,Ei,t) (12) 

The norm of Vj is a function of the energy of the particle given by the well 
known relation: 

Bk = f N 2 (13) 

Let us note that the distribution of density of mass pi, which can also be 
a sum of Dirac delta functions, depends on the angles 9, 0, on the position 
r of the point P in the space, on the energy E{ and on the time. There is 
however no dependence on the radial coordinate w. We will see below why 
it is not necessary to add the radial coordinate in the list of the arguments of 
Pi. Looking at Fig. El it is clear that the number of particles incident upon 
the surface dS from the specified direction is given by: 

dN i>eR {r, 9, cf>, E i: t) = ^| Vi |dS • e R (9, 4>)dt (14) 

In order to derive Eq. (jTlj) it has been used the fact that the total mass 
dM of the particles which traverse the surface dS in the interval of time dt is 
dM = Pi \vi\dS-e R dt. The number of such particles is obtained after dividing 
the total mass dM by the mass m ; of a single particle of type i. The depen- 
dence of dNi iBR (r, 9, <fi, Ei,t) on the kinetic energy E{ becomes explicit after 
eliminating the norm of the velocity |vj| from the right hand side of Eq. (jTll) 
using Eq. (JHJ). Finally, the scalar product dS ■ e R (8,cf>) gives the effective 
area which is hit by the particles incoming from the direction e R (9, <p). Since 
by definition of DDI the particle velocities are always perpendicular to dS, 
i. e. parallel to dS, see Eq. ([TO . we have that 



dN ie Jr,9,(P,E t ,t) = ^-WAdSdt = ^-J—dSdt (15) 

rrii rrii y 

Now we are in a position to understand why it is not needed to take into 
account the dependence on w of mass density. The reason is that just a small 
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Figure 5: This figure shows the geometrical setup for the definition of the 
differential directional intensity. The normal vector n to the infinitesimal 
surface dS at the point P coincides with the vector e R which gives the di- 
rection of the incoming particles. The particles that will be traversing the 
surface dS within the interval of time dt are those contained in the volume 
h ■ dS, where h = \\i\dt. 
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portion of space near the point P is considered, so that the radial coordinate 
is varying within the interval [0, h), where h = \Vi\dt. Clearly, the variation of 
Pi with respect to the radial coordinate is negligible within this infinitesimal 
interval. 

Let us remark that, in real measurements, the number of particles coming 
from a particular direction is usually very small, so that it is better to consider 
an entire set of directions, for instance those characterized by slightly different 
angles 9', 0' included within the range 

6<9'<9 + A9 (16) 
<P < <P' < ^ + A0 (17) 

where A9 and Acj) denote finite quantities and not infinitesimal ones. Clearly, 
the unit vectors eji(9' , <//) associated to these directions span a surface of area: 

A= / d9'sm(9') / d4' (18) 
Je J(j) 

on a sphere of unit radius, see Fig. 13 Always for experimental reasons, it 




Figure 6: Area A spanned on a sphere of unit radius by the unit vectors 
eji(9' ,<p') . The values of 9' and cf>' are defined in Eqs. (IT6l) and (TP?]) . 

will also be convenient to enlarge the set of possible particle energies to a 
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finite interval: 

E l <E' i <E i + AEi (19) 

Now we count the incoming particles in the neighborhood of the point P 
according to the following procedure. For each value of the energy E[ within 
the interval (fTH and for each of the directions falling within the range of 
angles of Eqs. (fT6HT7|) . we count the number of particles traversing a surface 
of fixed area dS and perpendicular to that direction in the sense of Eq. (fT0|) . 
Successively, summing over all the numbers obtained in this way, we get the 
final result dNi A 6 A< f >A E i { r ,0> ( j>, Ei,t)dSdt. In terms of mathematical equa- 
tions, this means that the quantity cL/V^A^A^A^r, </>, E^t) is computed 
by integrating the infinitesimal number of particles dN ijeR {r,9,(j),E^t) of 
Eq. JUJ) as follows: 

dN i , A 0,^ AEi {r,e^,E i ,t)= / / -^-jJ2Em'dSdt 

JEi Je J 4> (m^) 2 

(20) 

To write the above equation we have used Eq. (fl~3|) in order to express the 
speed |vj| as a function of the kinetic energy E^ and the fact that the in- 
finitesimal element of solid angle d£l! is given by: 

dQ' = sm(9')d9'd<f)' (21) 

The quantity (i/V^A^A^A^ 1 ', 0, E^t) describes the intensity of particles 
with energy in the interval (|19|) which arrive at the point P from all the 
directions spanning the area A of Eq. (|T£|) on a sphere of unit radius. 

To compute the DDI, it is now sufficient to take the limit in which AO 
and A</> become infinitesimally small. In that case: 

dN iAeA ^ AEi (r, 9, 0, E h t) ~ dNi(r, 9, 0, E it t) (22) 

with 

dNfa 9, 0, E u t) = —^J^EidttdEidSdt (23) 

K) 2 

It is easy to realize that the number of particles dN^r, 9, (j), E i: t) of Eq. 
coincides with the number of particles entering in the definition of DDI of 
Eq. (JHJ). Comparing these two equations, we find that 

I ddi (r, 9, <P, Et, t) = y^-r^Ei (24) 
[mi) ^ 
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Eq. (|2"Ij) provides a nice relation between the DDI and the mass density p». 
From Eq. (|23J) it turns out that the units in which the DDI is measured are 
cni~ 2 s -1 sr~ :L GeV~ 1 , where sr is a shorthand for steradiant, the unit of solid 
angles. In SI units one should replace centimeters by meters and GeV's by 
Joules. 

As a final remark, let us note that the DDI carries a lot of information 
about the intensity of particles incoming from different directions. We will 
also see that the DDI is related to the energy density of particles. On the 
other side, the DDI is an observable which is mainly related to the point in 
which it is measured. In particular, it is not possible to integrate the quantity 
Iddi in dS in order to extend its meaning to an arbitrary finite surface S. The 
reason is that, in the definition of the DDI the vector element of surface dS 
is constrained to satisfy Eq. (fTU|) . Of course, in a real measurement the 
infinitesimal element of surface dS is necessarily approximated by a finite 
surface AS, which may be for example the sensor of some particle detector. 
However, if one wishes to measure the number of particles traversing an 
arbitrary finite surface S, one should introduce the concept of flux. This will 
be done in Subsection IA.31 



A. 1.1 Quantities related to the differential directional intensity 

Starting from the differential directional intensity Iddi it is possible to con- 
struct several other quantities which are often encountered in the literature. 

Integral directional intensity: The integral directional intensity (IDI) I i( n 
is obtained by integrating the DDI over some finite interval of energy 
AE. — E- — E- ■ ■ 

Iidi(r,9, 0, E iimin ,E iimax ,t) = / I ddi (r,9,(p,E i ,t)dE i (25) 

where, of course, E i)min > in the non-relativistic case. Moreover, 
Ei t max £ [0,oo]. The units of I ic n are cm~ 2 s _1 sr _1 . 

Differential vertical intensity: In many experiments it is measured the 
intensity of particles arriving on the surface of the Earth from the 
vertical direction, i. e. the direction that from the ground points 
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toward the center of the earth 12 . For this reason, the DDI in the 
vertical direction, which is supposed here to coincide with the value of 
the angle 9 = of our spherical system of coordinates, has deserved a 
separate name and it is called the differential vertical intensity (DVI) 
Idvi- The Idvi is defined as follows: 

I dvi (r, E i} t) = I ddi (r, 9 = 0,(j), E h t) (26) 

The units of DVI are the same of the units of the DDI. 

Integral vertical intensity The integral vertical intensity (IVI) gives the 
number of particles coming from the vertical direction with respect to 
our coordinate system to, 9, <p and with energies comprised within the 
interval [E iiTnin , E i)max ] which traverse a unit surface in the unit of time: 

I im {r, t) = f El ' maX I ddi (v, 9 = 0,^, E h t)dEi (27) 

Integrated intensity: The integrated intensity (II) 13 In is defined as the 
integral of the DDI over all possible directions and energy values[39|: 

/*+oo p2tc /*7T 

I u (r,t)= dEi d<f) d9 sin 9 I ddi (r, 9, <P,Ei,t) (28) 
Jo Jo Jo 

The integrated intensity is measured in units cm~ 2 s~ 1 . Of course, if 
one integrates the DDI only over all possible directions the result is a 
quantity which may be called the differential integrated intensity. 

A. 2 Energy density of particles and the intensity in the 
relativistic case 

The specific directional density (SDD) u sdd (r, 9, (f>, Ei, t) is defined here as the 
kinetic energy of particles of type i per unit of volume, of energy and of solid 
angle. Sometimes the SDD is also called differential directional density. The 
quantity: 

dU = u sdd {r, 9, 4>, E h t)dVdVtdEi (29) 

12 Sometimes one considers the intensity of particles in the near vertical direction, where 
the values of the angle 9 between the trajectories of the particles and the gravity force 
spans over a finite interval, such as for instance 0.9 < sin0 < 1. 

13 Strictly speaking the name integrated integral intensity would be more correct. 
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represents the total kinetic energy carried by particles which are inside an 
element of volume dV and have velocities whose directions span a small 
element of solid angle dfl centered around the direction of the unit vector 
gr(0, <fi). The norms of these velocities are determined by the condition that 
the energy of the particles must be within the infinitesimal interval [Ei, Ei + 
dEi]. Clearly, the number of particles dN with the above characteristics 
which are inside the small volume dV at the time t is given by the total 
energy of the particles divided by the energy of each single particle: 

dE- 

dN = u sM (r,e,<f),E i ,t)dVdtt—^ (30) 

In the non-relativistic case we have seen that the relation between the kinetic 
energy Ei and the norm of the velocity |v»| is provided by Eq. (|T3|) . In the 
relativistic case, this equation must be substituted with the following one: 



Ei + rriiC 2 



Ef + 2m l c 2 E i (31) 



The relativistic and non-relativistic expressions of the energy density u sdd 
may be found in Ref.[40j. 

At this point, we wish to derive the DDI for particles which attain rela- 
tivistic speeds in terms of the energy density. We note to this purpose that 
the particles which will traverse the surface dS in the time dt while arriving 
from the direction perpendicular to dS, are contained in the infinitesimal 
volume: 

dV = \vi\dSdt (32) 

As a consequence, the number of particles dN t of Eq. © may be expressed 
in terms of the SDD as follows: 

dEi 

dNi = u sdd (r,0,(f),Ei,t)\vi\dSdtdn— - (33) 

Ei 

Comparing Eq. (J0J) with Eq. ijHHjl . we obtain a relation between the DDI and 
the SDD: 



hdi = u add (r, 6, 0, E u t)--^- ozry/fii + Irrii&Ei (34) 

where we have used Eq. ([3TJ) in order to write the speed |v$| as a function of 
the energy Ei. Eq. (|34| is the analogous of Eq. (|24| in the relativistic case 
with the mass density pi replaced by the energy density. 
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At this point, the derivation of the related quantities of the DDI, like for 
instance the integral directional intensity or the differential vertical intensity, 
proceeds as in the non-relativistic case of Subsection IA.11 

Finally, following an analogous calculation of the total energy density of 
electromagnetic radiation presented in Ref. [HE], it is possible to compute the 
total kinetic energy density u te d{r,t) of the particles per unit of volume: 

u ted {r,t)= dEi / dn—^= 1 1 I ddi {v,9,<P,E i: t) (35) 
J ™*i* J c^Ef + 2m i c 2 E i 

A. 3 Differential directional flux and related quantities 

The definition of the differential directional flux (DDF) $ddf is very similar to 
that of the DDL The difference is that, in the case of the DDF, the direction 
of the normal n to the element of surface dS does not need to coincide with 
the direction of the velocity of the incoming particles e R (9, (j)). More precisely, 
the DDF is defined in such a way that the quantity 

dN[ eR = <5> ddf dSdtdEidQ (36) 

represents the number of particles of a given kind traversing the infinitesimal 
surface element dS during the time dt within the element of solid angle dfl 
and within the energy interval [E i: Ei + dEi]. The superscript / has been 
added to remember that now a flux is being computed and not an intensity. 

To compute the DDF, let us imagine that we wish to measure it in a 
neighborhood of a point P. Since such measurements are usually performed 
on the ground, to fix the ideas we assume that the point P is very near 
(a few meters or less) to the surface of the Earth. The particle detector is 
approximated as a small and flat surface, which is centered around P. One 
side of the surface, that in which there are the sensors and which is thus 
able to detect the fluxes of incoming particles, is always directed toward the 
sky, while the opposite side is pointing toward the ground, see Fig. To 
express this all in mathematical terms, we introduce an infinitesimal element 
of surface dS and we choose a system of coordinates xp, yp, zp at the point P 
in such a way that dS lies in the horizontal plane zp = 0. In polar coordinates 
w, 6, (p this means that the direction of the unit vector n which is normal to 
dS is given by the angle 9 = 0. Furthermore, the orientation of n is such 
that it points downward, i. e. toward the ground. Since we are measuring 



25 



Sky 

/////////////////// 




Ground 



Figure 7: This figure shows the schematic experimental setup used to mea- 
sure the flux of particles of cosmic origin on the ground. The detector is 
represented as an infinitesimal element of surface dS. The active part of the 
detector is on the upper side of the surface, which points toward the sky. 



only particles which traverse the upper side of the surface dS in a downward 
sense, this implies that: 

n- e R {0,<j>) = cos 9 > (37) 

i. e. the normal vector n and the particle velocity |vj|e#(#,0) form an angle 
9, see Fig. El Clearly, Eq. f37l) is satisfied only in the interval < 9 < n2 
The volume of particles dV which will traverse dS coming from the direction 
e r (9, <fi) is shown in Fig. |E]and it is given by: dV = \vi\dt cos9. The number 
of those particle in the non-relativistic case is thus given by 

dK f e = ^^dSdt cos 9 (38) 
% ' BR rrii 

This is the analog of Eq. (fn| in the case of the flux. The factor cos 9 results 
after evaluating the scalar product n • e R (9,(j)) using Eq. (j3*TJl . As in the 
case of the intensity, in real measurements it is better to consider particles 
coming from different directions and carrying different energies, instead of 
fixing the attention to a particular direction and a particular energy. However 
we should remember that, in the computation of the DDF, the orientation of 
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Figure 8: This figure shows the geometrical setup for the definition of the 
differential directional flux. The normal vector n to the infinitesimal surface 
dS at the point P makes an angle 6 with the vector e R which gives the 
direction of the incoming particles. The particles of type % with velocity 
|vj|ejj(0, 4>) which will traverse the surface dS within the interval of time dt 
are those contained in the volume h ■ dS cos 9, where h = \vi\dt. 
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the surface dS remains fixed, what is changing is the direction of the incoming 
particles. At this point the procedure to obtain the differential directional 
flux is entirely similar to that used in deriving the explicit expression for the 
DDI of Eq. (125). The final result for the DDF is : 

$ ddf (r,6,(t),E i ,t) = -^J2E i cosd (39) 
(m^) 2 

The units of the DDF are the same as the units of the DDL 



A. 3.1 Observables related to the DDF 

In analogy with what has been done in the case of the intensities, one may 
construct other observables starting from the DDF. For example, after inte- 
grating the DDF with respect to the energy one obtains the integral direc- 
tional flux (IDF), while the differential vertical flux (DVF) corresponds 
to the value of the DDF in the case 9 = 0. Clearly, the DVF coincides 
with the DVI. The differential integrated flux $(r, E^ t), which can 
be found in standard textbooks, is defined in such a way that the quantity 
dN{ int = $(r, Ei,t)dSdt coincides with the number of particles of a given 
kind and of given energy Ei traversing in a downward sense 14 the element 
of surface dS: 

— 27T 

$(r,£i,t)= d6 J d0$ dd/ (r, #>,£,,£) sin 6> (40) 

The subscript int in dN( int means integrated and it is referred to the fact 
that, to derive $(r, Ei,t)dSdt, one needs an integration over dfl. This flux 
can be measured in units cm~ 2 s _1 GeV _1 . Starting from the expression of 
$(r, Ei,t) it is possible to compute the integral integrated flux, often called 
total flux: 

$(r,t) = f ^{r.E^dEi (41) 



The concept of flux is usually connected with vector fields. In the present 
case, the vector field is provided by the so-called differential directional in- 
tensity field 

Iddi = hdi^R (42) 

14 We recall that downward means here that the angle between the normal vector n and 
the vector e_R which determines the direction of the particle velocity must be within the 
range [0, £]. 
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where e R is the unit vector which defines the direction of the velocity of 
particles. In terms of the DDI field, the flux may be expressed as follows: 



according to the usual definition of flux. 

Contrarily to the DDI, the DDF may also be integrated with respect to 
the element of area dS. As we have anticipated before, this fact allows to 
define the flux of particles traversing an extended surface S. To compute the 
flux of particles in the case of an extended surface S, it will be convenient 
to parametrize this surface with the help of two parameters o\ and <r 2 , so 
that a point of S in the space will be denoted by the triplet of cartesian 
coordinates x(a±, a 2 ), y(cr 1 , cr 2 ) and z(a±, a 2 ) or, shortly, by the radius vector 
r(o"i,o"2)- For each point P of the surface, corresponding to a given value 
of the parameters (J\ and a 2 , we have seen that it is possible to compute 
the differential flux of particles <3>(r(o"i, 02), Ei, t) traversing a small element 
dS of S. The total flux &s(Ei,t) of particles of energy incoming upon 
S within the interval of time dt is obtained by integrating $(r(<7i, a 2 ), E { , i) 
with respect to dS, where dS will now depend on <i\ and a 2 : 
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